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We report an observation of a commensurate-incommensurate phase transition in the Dzyaloshinskii-Moriya 
spiral antiferromagnet Ba2CuGe2 07. The transition is induced by an external magnetic field applied along the 
c axis of the tetragonal structure, i. e., in the plane of spin rotation. Bulk magnetic measurements and neutron 
diffraction experiments show that the transition occurs in a critical field H c « 2.1 T. Experimental results for the 
period of the magnetic structure and magnetization as functions of magnetic field are in quantitative agreement 
with our exact analytical solution for Dzyaloshinskii's model of commensurate-incommensurate transitions in 
spiral magnets. 
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I. INTRODUCTION 

Among relatively unusual and exotic magnetic interac- 
tions in solids is the Dzyaloshinskii-Moriya (D-M) asym- 
metric exchange inter action. tMu Unlike the conventional 
Heisenberg exchange coupling, it is proportional to the 
vector product of interacting spins, and is permitted by 
symmetry only in non-centric crystal structures. It re- 
sults from relativistic spin-orbit cprrectionsS to the or- 
dinary superexchange mechanismu and therefore is usu- 
ally weak compared to antiferromagnetic symmetric ex- 
change. This is why only few materials in which D-M 
interactions play an important role have been f^wid so 
far 
MnS 



best known examples are cubic FeGeO'El and 
where D-M terms in the Hamiltonian cause an 
instability of ferromagnetic order towards the formation 
of an incommensurate spiral structure. 

This paper deals with the properties of Ba2CuGe 2 07, a 
newly investigated system in which D-M coupling plays 
a key role. Originally the magnetic properties of this 
quasi 2-dimensional (2D) antiferromagnet (AF) were in- 
vestigated as part of the ongoing search for new singlet 
ground state compounds, triggered by the discovery of 
a spin-Peierls transition andj-oLhcr extraordinary mag- 
netic properties of CuGe03.E£lEil While some .new ma- 
terials related to CuGeC>3, e.g., CaCuGe206r-3il3 and 
BaCuSi2 06,til have dimerized ground states and energy 
gaps in their spin excitation spectra, Ba 2 CuGe2 07 un- 
dergoes a transition to a magnetically ordered phase be- 
low Tjy ~ 3.2 K_and the magnetic excitations are gap- 
less spin waves.tJ Even though the magnetic properties 
of Ba2CuGe207 can be adequately described in terms of 
classical spins, they are rather intriguing: the magnetic 
structure is an incommensurate spin-spiral. By perform- 
ing detailed measurements of the spin-wave dispersion 



we have previously demonstrated that this spiral orderiiag 
may not be caused by competing exchange interactions Jij 
Since Ba2CuGe207 has a non-centric crystal structure, 
we have suggested that the incommensurate magnetic 
phase is a result of Dzyaloshinskii-Moriya interactions. 

The crystal structure of Ba 2 CuGe207, as well as 
the spiral spin arrangement in the magnetically ordered, 
phase, were discussed in detail in our previous work,li3 
and only the essential features are reviewed here. The 
magnetic (S — 1/2) Cu 2+ ions are arranged on a square 
lattice in the (a, b) plane of the tetragonal structure 
[space group P42im (No. 113), lattice constants a — 
8.466 A, c = 5.445 A]. The in-plane interactions between 
spins are established through GeC>4 tetrahedra that, to- 
gether with the Cu 2+ ions, form distinct Cu-Ge-0 lay- 
ers. These layers are well isolated by interstitial non- 
magnetic planes of Ba 2+ ions. Only nearest-neighbor 
in-plane antiferromagnetic exchange interactions are im- 
portant (J s» 0.48 meV). Interplane coupling is ferromag- 
netic and substantially weaker (Jj_ w 0.013 meV). In the 
ordered phase the spins lie in the (1, 1, 0) plane and the 
propagation vector for the spiral is (1 + £, £, 0) where 
C = 0.027 [Fig. 0(a)]. The magnetic structure is a distor- 
tion of a Neel spin arrangement: a translation along the 
(1,1,0) direction (hereafter referred to as the x axis) in- 
duces a rotation of the spins by an angle a = =£■ « 9.7° in 

the (1, 1,0) plane (relative to an exact antiparallel align- 
ment). Along the (1, 1,0) direction (y-axis) the spins are 
perfectly antiparallel. Nearest-neighbor spins from adja- 
cent Cu-planes are aligned parallel to each other. 

The mechanism by which D-M interactions can sta- 
bilize the spiral structure in Ba2CuGe207 is illustrated 
in Fig. |l|(b). The D-M energy for two interacting spins 
Si and S2 may be written as (Si x S2) • D^ 1,2 \ where 
D^ 1 ' 2 ) is the so-called Dzyaloshinskii vector attributed to 
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the oriented bond between the two spins and "x" de- 
notes a vector product. For spins lying in the (1,1,0) 
plane the only relevant component of the vector D is 
D y . The symmetry of the structure is such that for two 
subsequent Cu-pairs (1,2) and (2,3) along the (1,1,0) 
direction D y 1,2] = D y 2,3] [Fig. 0(b)]. A finite value for 
D energetically favors a 90° angle between subsequent 
spins. Since AF exchange favors an angle of 180°, the to- 
tal energy is minimized at some intermediate angle that 
is defined by the ratio of exchange and D-M interactions 
strengths. 

How can the mechanism described above be verified 
experimentally? If one could force the spins into the 
(0, 0, 1) plane for example, D y would become inactive, 
due to the nature of the vector product, and the incom- 
mensurate spiral would disappear. The z-component of 
D, i.e., the (0, 0, l)-projection, on the other hand, would 
become relevant. As shown in Fig. 0(b), for the D Zl 
unlike for D y , symmetry dictates a change of sign from 

one Cu-Cu bond to the next: = —D^' 3 \ As a 

result, D-M interactions would tend to distort the Neel 
spin arrangement towards a weak-ferromagnet (canted) 
structure, that is commensurate with the crystal lat- 
tice. In practice forcing spins into the (a, b) plane can 
be achieved by applying an external magnetic field along 
the c-axis of the crystal. Indeed, in structure with zero 
net magnetic moment spins tend to align perpendicular 
to the external field, which in the simplest case leads 
to spin-flop transitions in conventional antiferromagnets. 
In other words, if the proposed model for D-M interac- 
tions in Ba2CuGe207 is correct, we expect a field-induced 
commensurate-incommensurate magnetic transition (CI) 
in this material. 

In the present paper we report an experimental ob- 
servation of such a phase transition in Ba2CuGe207 by 
means of bulk magnetization and neutron diffraction 
measurements. We find that the phase transition is of 
rather unusual character and is the first "clean" realiza- 
tion of Dzyaloshinskii's model for CI transitions in spiral 
magnets induced-by a magnetic field applied in the plane 
of spin rotationll3 We also present an exact continuous- 
limit solution to Dzyaloshinskii's model. For the partic- 
ular case of Ba2CuGe207 theoretical predictions are in 
quantitative agreement with experimental data. A shoxt 
preliminary report on this work is published elsewhere.EZl 

II. EXPERIMENTAL PROCEDURES 

Transparent, slightly yellowish single-crystal samples 
were prepared using the floating-zone method. Magne- 
tization measurements were performed with a conven- 
tional DC-squid magnetometer in the temperature range 
2-300 K. Two sets of neutron scattering experiments were 
carried out on the H4M (thermal beam) and H9 (cold 
beam) 3-axis spectrometers at the High Flux Beam Re- 
actor (HFBR) at Brookhaven National Laboratory on an 



irregularly shaped sa 4 x 4 x 4 mm 3 single-crystal sam- 
ple with a mosaic spread of » 25'. In the first set of 
measurements on H4M (experiment I) the H-T phase 
diagram was determined using neutrons of incident en- 
ergy Ei = 14.7 meV with a (20' - 40' - 20' - 40') col- 
limation setup and two PG filters. The second set of 
measurements was done on H9 (experiment II) and the 
field-dependence of the magnetic propagation vector was 
studied using (60' - 40' - 60' - sample - 15' - 80') colli- 
mations, an Ei = 4.6 meV neutron beam and a Be-filter 
in front of the sample. In both experiments the use of a 
pumped 4 He cryomagnet allowed us to work in the tem- 
perature range 1.3-5 K and magnetic fields up to 6.5 T. 
The sample was always mounted with the (h, k, 0) zone 
in the scattering (horizontal) plane and the field was ap- 
plied along the vertical direction (c-axis of the crystal). 

The sample we used in preliminary experiments shat- 
tered when it was cooled down in the cryomagnet for 
the second time. We therefore tried to mount the new 
crystal as strain-free as possible. The new sample was 
wrapped in aluminum foil that was attached to a thin 
Al-plate. This technique has one serious drawback. As 
will be explained below, a good alignment of the crystal- 
lographic c axis with the magnetic field is crucial. Un- 
fortunately, the described mounting procedure does not 
allow to maintain an alignment of better than w 1°. Be- 
fore the sample is put into the cryostat, an almost perfect 
alignment is achieved. It is upon cooling to base tem- 
perature that the undesirable misalignment occurs. In 
experiment I the (1,-1,0) and (1,1,0) crystallographic 
directions formed angles of 1.5° and 1° with the horizon- 
tal plane, respectively, as measured in situ at low tem- 
perature. In experiment II these angles were 1° and less 
than 0.3°, respectively. 

III. EXPERIMENTAL RESULTS 

A. Magnetization 

The first evidence of a field-induced magnetic phase 
transition in Ba2CuGe207 was found in bulk magnetic 
measurements. The longitudinal magnetization data col- 
lected at T = 2 K is plotted against magnetic field H 
applied along the a or c crystallographic axes in the in- 
sert in Fig. 0. For H||a no anomalies are observed. In 
contrast, when the field is applied along c, M(H) has a 
broad step-like feature around H = 2 T. The anomaly 
is best seen in the plot x(H) that was obtained by nu- 
merically differentiating the experimental magnetization 
curve (Fig. 0, main panel). 

B. Neutron diffraction 

The phase transition is best observed in neutron 
diffraction experiments. Figure (experiment II) shows 
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some elastic scans along the (1 + CiCjO) direction mea- 
sured in Ba2CuGe207 at T — 1.4 K. The different scans 
correspond to different values of magnetic field applied 
along the c-axis. In zero field rfFig- the magnetic 

peak is positioned at £ = 0.027.E-3 The incommensurabil- 
ity parameter £ decreases gradually with increasing H. 
Around H = 1.8 T an additional peak appears at the 
Neel point (1, 0, 0). At the same time the intensity of the 
satellite at (1 + £, £, 0) starts to decrease rapidly and the 
rate at which £ changes with field is increased [Fig. §(b)]. 
On further increasing the field the central peak gains in- 
tensity and the satellite eventually vanishes at some crit- 
ical field H c S3 2.3 T [Fig. |(c)]. At higher fields only the 
Neel (commensurate) peak is seen, all the way up to the 
highest field available experimentally. The intensities of 
the (1 + OCiO) satellite and the (1,0,0) peak are plot- 
ted against temperature in Fig. ||(d). Measured C(T) is 
shown in solid circles in Fig. ||. 

C. Domains and sample alignment 

The field-dependent behavior is extremely sensitive to 
the alignment of the crystallographic c axis with the 
applied field. The limitations of the sample-mounting 
technique employed, together with technical impossibil- 
ity of adjusting the alignment in situ at low temperature, 
are therefore a serious experimental complication. This 
problem is closely linked to the issue of magnetic do- 
mains. In the tetragonal symmetry two spiral structures 
with propagation vectors (1 + £, 0) an d (1 + d — Cj 0) 
are allowed. Zero-field cooling produces domains of equal 
volume. If the sample is cooled in an H = 1 T field, that 
is removed only at low temperature, only one domain is 
present. Since we are able to exactly measure (but not 
adjust!) the orientation of the sample in situ, we could 
determine that it is the domain for which the plane of the 
spiral forms a larger angle with the applied field that is 
stabilized. This behavior can easily be explained. If the 
plane in which the spins are confined is not parallel to the 
external field, the system can gain some Zeeman energy 
without sacrificing much D-M or exchange energy: all 
spins may tilt slightly in the direction of the field compo- 
nent normal to the spin plane, producing a cone structure 
with practically no change in the angles between neigh- 
boring spins. The domain for which the spin plane forms 
a larger angle with the applied field is thus energetically 
more favorable. Note that in the commensurate (high- 
field) structure there are also two possible domain types, 
with spins in the (1, 1, 0) or (1, 1, 0) planes. 

The data shown in Fig. were collected on a field- 
cooled sample and the satellites correspond to the more 
misaligned (1 + £,00) domain (domain A, S3 1° tilt). 
No satellites were observed in the (1 + £, — £,0) domain 
(domain B, < 0.3° tilt) in any fields in field-cooling ex- 
periments. When the measurements at T = 1.4 K were 
done on a zero field cooled sample, we observed that scat- 



tering intensities originating from both domains remain 
comparable in fields up to S3 1.7 T. At higher fields the 
Neel peak appears and domain B is rapidly destroyed. 
Applying a sufficiently large field thus brings us back to 
the situation where the entire crystal is a single magnetic 
domain. 



D. Phase diagram 

The experimental H — T phase diagram for 
Ba2CuGe2 07 is shown in Fig. ||. The I-O-II line cor- 
responds to the appearance of the magnetic reflection at 
the commensurate (1,0,0) position. The III-O-IV line 
shows where the satellite peaks disappear. Both the O-II 
and the O-IV lines were measured in a field-cooled sam- 
ple and the misalignment of the corresponding domain 
was S3 1.5° (experiment I). 

In the high-field phase the structure is commensurate, 
with a propagation vector (1,0,0). This may be the 
signature of Neel order or (see introduction) a canted 
"weak ferromagnetic" structure. In the low-field region 
the structure is incommensurate with a propagation vec- 
tor (1 + £, £, 0). Only one domain is present since, as 
mentioned, field-cooling was used. Note that except for 
the incommensurability of the low-field phase, the phase 
diagram strongly resembles that of an easy-axis Heisen- 
berg antifcrromagnet in external field applied along the 
easy axis. 

In the II-O-IV region the Neel peak coexists with the 
incommensurate satellite. It would seem that in the mis- 
aligned domain the field-induced transition is of first or- 
der and that the II-O-IV region corresponds to a mixed- 
phase state. Indeed, at T = 1.4 K we have observed some 
field-hysteresis: on sweeping the external field down from 
H = 4 T the Neel component vanishes and the satellite 
appears at slightly lower fields, than in the case when the 
H is gradually increased from zero. We believe that in 
a perfectly aligned sample the transition must be of 2nd 
order, in which case the II-O-IV region is replaced by a 
single line. 

E. Measurements of C(-ff) 

Accurate measurements of the field-dependence of £ 
were performed in experiment II where one domain was 
tilted by S3 1° (domain A) and the other almost perfectly 
aligned (domain B). Measurements for domain A could 
easily be done at base temperature (1.4 K) in a field- 
cooled sample (Fig. ||, solid circles). Additional data for 
domain A were collected at T = 2.4 K. At this tem- 
perature £(-ff) was found to follow the same curve as at 
T = 1.4 K to within experimental error. 

For domain B experiments in the interesting field re- 
gion, where £ starts to decrease rapidly with increasing 
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H, could not be performed at low temperature: the do- 
main itself is destroyed in fields higher than k 1.7 T (see 
above). Nevertheless, ((H) f° r domain B could be mea- 
sured in a zero-field-cooled sample at T = 2.4 K, i.e., just 
at the temperature of the critical point O in the phase 
diagram. At this temperature the B-domain survives up 
to the critical field and the Neel component does not ap- 
pear before the (1 + £, — £, 0) peak vanishes. Apparently 
even at 2.4 K domain wall pinning is sufficiently strong 
and kinetics sufficiently slow to make the energetically 
less favorable (better aligned) domain metastable even 
in high fields. This does not happen at low tempera- 
ture: the commensurate structure emerges in a lst-order 
transition and presumably the thermodynamic force that 
destroys the better aligned domain is larger. Some typi- 
cal scans for domain B at T — 2.4 K are shown in Fig. ||. 
The dashed line shows experimental Q-resolution. The 
(1 +Ci — (, 0) satellite is resolution-limited at H < 1.95 T 
[Fig. |q(a)] and starts to gradually broaden at higher fields 
[Fig. Mb, c)]. The experimental (,{H) is plotted in open 
circles in Fig. || We clearly see the sensitivity of the sys- 
tem to even a minor tilt of the applied field relative to 
the c axis. While at low fields ((H) for domains A and 
B coincide, there is a substantial discrepancy near the 
transition point. 

F. Critical indexes 

In our previous work we have reported measure- 
ments of the order-parameter critical exponent /3 = 
0.145(0.005) for the zero-field transition from paramag- 
netic (PM) to spiral states. Having discovered the com- 
mensurate phase, we have in addition determined /3 for 
the PM — > commensurate transition at H = 3.66T. The 
intensity of the (1,0,0) magnetic reflection is plotted 
against (Tjy — T)/Tjv in a log-log plot in Fig. |?] (open 
circles). Fitting the data with a power-law dependence 
we obtain Ph=4T = 0.185(0.005). The previously mea- 
sured temperature dependence of the (1.0273,0.0273,0) 
Bragg intensity at H = is shown in solid symbols in 
Fig. u[ Although the experimentally determined critical 
indexes are rather close to each other, in our data we 
clearly see that (3 = 0.145 is incompatible with the 4 T 
measurement (Fig. ^, dashed line). 

IV. THEORY 

In the introduction section we have qualitatively shown 
how D-M interactions stabilize the magnetic spiral in 
Ba2CuGe2C>7 and why applying a field along the c axis 
of the crystal should result in a commensurate structure. 
We shall now develop a quantitative description of a D-M 
spin-spiral in an external magnetic field. The general ex- 
pression for the free energy of such system was obtained 
by Dzyaloshinskii Jij who also predicted the CI transition 



and obtained analytical expressions for the critical prop- 
erties [at (H c — H)/H « 1]. In the present work we 
extend Dzyaloshinskii's approach and derive analytical 
expressions for the period and the uniform magnetiza- 
tion of the spiral valid throughout the phase diagram. 
We start by considering a simple, yet illustrative, case of 
a classical spin chain at T = 0. 

A. Classical spin chain with D-M interactions. 

Let us consider the 1-D case: a uniform chain of clas- 
sical spins with isotropic (Heisenberg) antiferromagnetic 
exchange. In addition, we include a Dzyaloshinskii en- 
ergy term with the D vector pointing along the arbitrary 
chosen y-axis. The system is then characterized by the 
following Hamiltonian (energy functional): 

W = 2J^S„S„ +1 +D^(S 

n X Sn+l)y (1) 

n n 

It is straightforward to show that to take full advantage 
of the Dzyaloshinskii term, the spins have to be confined 
to the (x, z) plane. The total interaction energy of a 
pair of spins is given by 8 — 2JS 2 cos <f> + DS 2 sin </> = 
5 2 V4J 2 + D 2 cos(<fr — a), where (f> is the angle between 
the spins, and a — arctan_D/2J <C 1, provided J ^S> 
D. The classical ground state is therefore an AF spiral, 
where the angle between subsequent spins is equal to 7r + 
a. 

Physically interesting behavior occurs when an ex- 
ternal magnetic field is applied perpendicular to the 
Dzyaloshinskii vector D, e.g., along z-axis. In this 
case one has to take into account the Zeeman energy 
—HgiiB^2 n S z , where g is the gyromagnetic ratio and 
/is is the Bohr magneton. Denoting by (f> n the angle the 
n-th spin S n forms with the z-axis, we can rewrite the 
energy functional as: 

U = ^ [ J cos(^„ + i -<j> n — a)-H cos <j> n ] , (2) 

n 

Where J = S 2 ^AJ 2 + D 2 and H = Hgfi B S. 

In Eq. (|^) we clearly see the competing terms that drive 
the CI transition. The J-term favors an incommensu- 
rate (spiral) structure, where the phases <p n of subse- 
quent spins differ by tt + a. The effect of the mag- 
netic field is more subtle. For H 3> J the spins tend 
to align parallel to the field. For H -c J however, the 
system still can gain some Zeeman energy without sacri- 
ficing much of its exchange energy. This is achieved in a 
spin-flop state, where all spins are roughly perpendicular 
to the magnetic field and slightly tilted in its direction. 
The external field thus favors a commensurate spin-flop 
structure, but it can only be realized by sacrificing some 
Dzyaloshinskii energy. Indeed, the energy difference be- 
tween the spin-flop configuration and the spiral state is 
J(l — cos a) — H 2 /8J. The spin flop state becomes en- 
ergetically favorable in fields H > H c ~ 2a, J . This is 
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the simplified physical picture already discussed in the 
introduction section. The above crude arguments does 
not take into account the distortion of the spiral by the 
applied field. The exact result for our model, H c = irJa 
is derived below. 



B. Exact ground state in the continuous limit 

It is convenient to change angle variables so that 

(j) n =7771 + 6 n (3) 

The energy functional in this notation becomes 
H = [-Jcos{6 n+ i -0 n -a)- (-l)"#cos0 n ] . (4) 

n 

The advantage of the new set of variables is that assum- 
ing a = 2ir( << 1 and H « J, the phase difference 
&n+i — n — a « 1, so we can safely replace the 1st term 
in Eq.(|) with -J + J(0 n+1 - n - a) 2 /2. 

The ground state satisfies the extremal conditions 
dTL/dOn = 0, which gives us the following set of equa- 
tions: 



H c = iraJ = 2w 2 (J 



'n+l 



29„ 



^(-l) n sin( 
J 



(5) 



These equations can be solved in the continuous limit. 
We shall look for a solution in the form n — Q(n) + 
(— l) n, d(n), where both 8(n) and #(n) are functions that 
only slowly change on the scale of a single lattice spac- 
ing. From Eq. (^) we obtain the extremal condition that 
should be satisfied in an energy minimum: 



d 2 e 

dn 2 



H 2 
— — sin 29 

8J 2 

■3 = -^sin9 
4J 



(6) 
(7) 



The first of these equations is the famous sine-Gordon 
equation, which allows for exact soliton solutions. The 
trivial "no-soliton" solution to Eqs. (|J|7]) is 6(n) = 7r/2, 
and i)(n) — —H/4J. This is precisely the spin-flop phase, 
realized in strong enough magnetic field. Less trivial is 
the one-soliton solution: 

6(n) = 2 arctan[exp(n_ff/2J)] - tt/2 (8) 
tf(ra) = -H/4Jt&nh(nH/2J) 

which has boundary conditions 0(— oo) = — ir/2 and 
9(+oo)=tt/2. 

For a — the energy of the one-soliton state will be 
always higher than that of the soliton vacuum (spin-flop 
state). In the general case, substituting (||) into the en- 
ergy functional (^J) , we find the energy difference between 
the one-soliton and soliton- free state: £\ — £q = H — Jira. 
We see that for H < H r , where 



H c = 4ir 2 (- 



JS 
9Hb 



(9) 



the energy of a single soliton is negative and solitons 
spontaneously "condense". This process will eventually 
saturate, due to the mutual repulsion of solitons, and the 
ground state will be a periodic "soliton lattice" . 

At the non-negligible soliton density the interaction 
starts to change the shape of individual solitons and has 
to be taken into account explicitly. Fortunately the gen- 
eral solution of the sine-Gordon equation is known: 



O(n) 



f3dx 



sin x 



nH 
~2J 



H 



0(n) = 



0(n) = amfn — —,0) 
2J(3 

— = sin Q(n) = =sn(n — —, 0) 

4J 4J 2J0 



(10) 

(11) 
(12) 



where om(i, 0) and sn(x, 0) are Jacobi elliptic functions 
of modulus 0. Substituting this solution into the energy 
functional (^|) we obtain the interaction energy per spin: 



£ 



16J 



,2 r 



2E{0) 



2H r 



2 2 K{0) 0HK(0) 



(13) 



where K and E are the complete elliptic integrals of the 
first and second kind, respectively. 

Equation (|l3|), as well as the expression for H c [Eq. (||)] 
were derived in Ref. |l^. In the present work we go one 
step further and determine which of the solutions given 
by Eqs. ( |ll| , |l2"| ) indeed correspond to the ground state 
at any particular H < H c . Unlike that in Ref. [l6], our 
treatment is not limited to the narrow critical region 
(H c — H)/H c « 1. Taking the partial derivative of 
the energy functional ( |l3| ) with respect to 0, andpSjibsti- 
tuting identities for derivatives of elliptic integralsc3 after 
some algebra we get: 



H 



P 



£ 



E{0) 
(9VbH) 2 1 
16J 2 



(14) 
(15) 



Equation ( |14[ ) is to be solved with respect to for any 
given H / H c . Q{H) can then be expressed in terms of 0: 



H TT 2 



C(H) 

C(0) H c 40K(0) 4E(0)K(0) 



(16) 



Equations ([14|) and (^6|) provide us with ((H) in a para- 
metric form. The ratio £(7J)/£(0) is plotted against 
H/H c in Fig. |[ The insert schematically shows the spin 
structure for H = 0, < H < H c , and H > H c . 

At low density soliton repulsion is exponentially small 
and the system close to the transition point is extremely 
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sensitive to any perturbations of the original Hamilto- 
nian, that will tend to pin down the soliton lattice in 
one of the metastable configurations. As H — > H c — 0, 
the solitan density decreases as l/ \ \n(H c — H)\, i.e very 
rapidly.La The transition is thus almost first order [the 
critical exponent S = in ((H) ~ (H c — H) s }. 

The soliton lattice at H ^ 0, unlike the pure sinu- 
soidal spiral at H = 0, has higher Fourier harmonics 
at C2n+i = ±(2n + 1)((H). The knowledge of the ex- 
act ground state enables us to calculate these harmonics 
analytically. For the Bragg intensity (square of the am- 
plitude) of the strongest 3rd harmonic, for example, one 
gets 



I 3 (H) q 2 (l- q 2 + q ±) 
h(H) (l +(7 2 + (Z 4)2 

I 3 (H) 2^(^/(1 -q 3 ) 2 + q 3 /(l + q 3 ) 2 ) 



(3 2 K(fif 



(17) 
(18) 



, where q = exp(-?r -^—^ ) 



This ratio of intensities is plotted against ((H) /((0) in 
Fig. |[ Except very close to H c all harmonics are much 
weaker than the first one. 

As we derived a parametric expression for ((H), so we 
can also obtain an exact parametric form for the uni- 
form magnetization M = —d£/dH. From Eqs. (H^H); 
and the identities for derivatives of elliptic integrak-Ej one 
gets: 



M 



(wb) 2 H 1 
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(19) 



The magnetization curve is continuous at the critical 
field. On the other hand, its derivative with respect to 
H (the field dependent susceptibility) diverges when H c 
is approached from below. 



C. Generalization to finite temperatures and higher 
dimensions 

The results obtained in the previous chapter can be 
easily generalized to the cases of higher system dimen- 
sionality and finite temperatures. We follow the approach 
developed by Dzyaloshinskii in Ref. |l6| and introduce the 
expression for the free energy of a general almost anti- 
ferromagnetic spiral in an external magnetic field. We 
generalize the results derived for the 1-D chain in the 
previous section to obtain analytical expressions for the 
period and the magnetization of the structure valid not 
only in the critical region, but troughout the phase dia- 
gram. 

Since the period of the spiral is long, nearest neighbor 
magnetic moments are almost opposite to each other. 
This enables us to approximately describe the structure 
in terms of position dependent staggered magnetization 



L(r). To take full advantage of the D-M interaction, L(r) 
must be perpendicular to the D-M vector, and therefore 
be confined to the x — z plane. We further assume that 
the magnitude |L(r)| does not depend on r. This sim- 
plification is justified since the Zeeman energy, which at 
H c is of order of D-M interaction energy, is much smaller 
than the energy of isotropic exchange. The local mag- 
netic structure is thus fully defined by the angle 6(r), 
that L(r) forms with the z-axis. The total "elastic" free 
energy per single Cu plane, associated with spin-spin in- 
teractions, can then be written as: 
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Equation (|20|) is given in the form suitable for describ- 
ing a single Cu plane in Ba2CuGe207. p s (T) is the so- 
called spin stiffness at given temperature, A = A a j, is the 
in-plane nearest neighbor Cu-Cu distance, and 7 is the 
spin stiffness anisotropy factor, defined as 7(A c /A Q f,) 2 = 
|«7 c j/|«7 a f,| » 1/37 in Ba 2 CuGe 2 07. Coordinates x, y, and 
z run along (1,1,0), (1,1,0), and (0,0,1) directions, respec- 
tively. The spuj-stiffness at zero temperature is given by 
p s (0) = 2JS 2 tB As T approaches Tjv the spin stiffness 
decreases as |L(r)| 2 and vanishes precisely at Tjy. 

In taking into account the external magnetic field, 
it is convenient to introduce the quantities x\\ (T) and 
X±(T), the magnetic susceptibilities of the system with 
respect to a field that rotates along with the spiral 
structure and is always parallel or perpendicular to the 
local staggered magnetization, respectively. As in a 
usual antiferromagnet, these susceptibilities do not di- 
verge at Tjv. In the paramagnetic phase (at T > Tjy), 
X||(T) = x±(T). At zero temperature the classical result 
is X±(0) = (.9/i B )7(8d/A 2 )@ and X \\(0) = 0. 

The full expression for the free energy density in the 
presence of magnetic field H applied along z-direction is: 
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At T — and d — 1 Eq. ( pT| ) coincides with the expression 
for the energy of a spin chain derived in the previous 
section. 

The equilibrium configuration of 6(r) minimizes the 
free energy and therefore satisfies a generalized form of 
Eq. (0): ' 



dx 2 



(X±~X\\)H 2 



sin v cos f = 



~ sin 26, (22) 
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where V — (p s /H 2 (x_L — Xll)) 1 ^ 2 - The critical field is now 
given by 

For a 4-nn AF like Ba2CuGe 2 0y at T — we can use 
X± = (gfi B ) 2 /l6JA 2 7 X|| = 0, and p.. = J ab S 2 . Substi- 
tuting these expressions into Eq. (B3h we obtain 

jq 

H c = 4V2(ir 2 — , (24) 

a factor of \/2, compared to the 1-D case [Eq. (g)]. 
This factor is due to the fact that in 2 dimensions 
each spin is antiferromagnetically coupled to four, rather 
that two nearest neighbors, and the overall structure is 
stiffer. Unlike antiferromagnetic in-plane interactions, 
ferromagnetic coupling between adjacent Cu-planes in 
Ba2CuGe207 does not change x_l or the expression for 
the critical field. 

The temperature dependence of H c can be easily un- 
derstood within the mean-field (MF) approximation. In 
this framework x_L is T-independent, x\\ decreases with 
Tn—T and x± = X\\ at T = TV. The effective strength of 
exchange coupling, represented by JS 2 in Eq. ( p3| ) goes 
as the square of the order parameter, i.e., as T/v — T (the 
MF order-parameter critical exponent (3 = 0.5). Substi- 
tuting these values in Eq.(|2^), we find that within the MF 
approximation H c is independent of temperature. This 
result is the same as for the spin-flop field in a conven- 
tional easy-axis Heisenberg antiferromagnet. 

Finally, we can generalize the expression ([l9]) for the 
magnetization. With new parameters in the sine-Gordon 
equation we get: 

for H < H c . For H > H c one has M = x±H. The field 
dependent susceptibility x{H) = dM/dH is given by: 

(m = | (XJ- -X||) / E{f3f 2E{0)\ 

X{ ' X[l [3 2 \ {l-f3 2 )K{l3f K{[3) ) 

(26) 

At H = this expression gives x(0) = (x_L + X||)/2: the 
structure is a uniform spiral, which effectively averages 
out the susceptibility for all directions in the spin plane. 
The-ausceptibility diverges at H c as 1/(H C — H) hi 2 (H c — 
H) x3 Above H c it has a constant value equal to x_i_ ■ 

V. DISCUSSION 

A. D-M interactions in Ba2CuGe2C>7 

We have shown that the spiral spin arrangement in 
Ba2CuGe 2 07 is due to the in-plane component D y of 



the Dzyaloshinskii vector D for nearest-neighbor Cu 2+ 
ions. It is rather difficult to experimentally determine 
whether or not the out-of-plane component D z is also 
active. In the spin-flop phase the predicted canted weak- 
ferromagnetic structure differs from a Neel state only in 
that it gives rise to additional magnetic peaks coinci- 
dent with nuclear Bragg reflections. The latter are much 
stronger than any magnetic scattering intensities, and 
make measurements of the ferromagnetic component all 
but impossible. In the spiral phase however, if D z =/= 0, 
satellites of type (h + (, k + £, I) should be observable 
around ferromagnetic zone-centers, in addition to the 
principal magnetic peaks positioned around the AF zone- 
centers. In preliminary experiments at T — 1.5 K we 
have indeed observed extremely weak elastic features at 
reciprocal-space positions (2 + £,£,0) and (1 + C, 1 + C; 0). 
So far we have not investigated the possibility of these 
features being artifacts due to nuclear-magnetic double 
scattering.Ej We plan to resolve this uncertainly in fu- 
ture experiments. 

B. The commensurate-incommensurate transition 

Studies of CI phase transitions have a long history, daU 
ing back to pioneering works ofjErenkel and KontorovaEil 
and Frank and van der MerweJ23. Since then CI transi- 
tions were discovered and studied experimentally in a 
number of such seemingly unrelated systems— as noble 
gas monolayers adsorbed on graphite jSurfacecJ, charge 
density Bzave materialscJ, ferroelectricscJ and rare-earth 
magnetsEj (For comprehensive reviews see for example 
Ref.E3). As a rule, CI transitions result from a com- 
petition between two distinct terms in the Hamiltonian 
that have different "built-in" spatial periodicities and are 
often referred to as potential and elastic energy, respec- 
tively. The potential energy by definition favors a struc- 
ture commensurate with the crystal lattice. Such is the 
interaction between gas atoms and the graphite matrix 
in intercalated and adsorbed systems. The elastic term is 
intrinsic to the system where the transition occurs, and 
has a different "natural" built-in period. For adatoms on 
graphite this term represents their mutual interaction. 
In our case of a Heisenberg AF with D-M interactions 
it is the H cos (f> n term in Eq. (0) that plays the role of 
an effective potential, forcing the spins in the plane, and 
thus favoring a commensurate structure. The competing 
elastic term is Jcos(0„+i — 4> n — a), and the "natural" 
periodicity is set by the angle a. 

In many known realizations of CI transitions, such as 
adsorbed gas monolayers, it is the period set by the elas- 
tic term that can be varied in an experiment to drive 
the transition, whereas both the strength and the pe- 
riod of the potential remain constant. In other sys- 
tems, among them rare-earth magnets, both the elastic 
term (exchange coupling between spins) and the poten- 
tial (magnetic anisotropy) can be changed, but only in- 
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directly, by varying the temperature. In both cases one 
typically observes a "devil's staircase" phase diagram (for 
a review see for example Ref. |27j): the incommensurate 
structure tends to lock onto rational fractions of the pe- 
riod of the potential. Instead of a continuous CI tran- 
sitions one gets a series of commensurate-commensurate 
transitions between different lock-in states. 

The interest of Dzyaloshinskii's model for CI transi- 
tions is that it is driven by a changing strength of the 
potential alone, with both built-in periods remaining con- 
stant. The experimentalist has a convenient handle on 
the potential term that he can vary by simply adjust- 
ing the external field. Since the "built-in" periods do 
not change, the transition is continuous with no "devil's 
staircase" behavior. 

The other advantage of the present realization of CI 
transitions lies in the fact that the potential energy has a 
pure sinusoidal form as a function of the angle 6(x). The 
model can be exactly solved and a quantitative compar- 
ison of theory and experiment are possible for the entire 
phase diagram. In most other CI systems one knows only 
that the potential energy is a periodic function with a 
given period, determined by the underlying lattice, while 
its exact functional form remains undetermined. Quanti- 
tative comparison of experiment and theory in this case 
is restricted to the narrow critical region close to the 
transition point. 

It is important to note here that while field-induced 
CI transitions have been previously observed in a 
number pf. magnetic insulators with incommensurate 
structurecaSEJfl for various reasons none of them is well 
described by Dzyaloshinskii's model. In FeGe (Ref. |]) 
the problem is that in the high-symmetry cubic struc- 
ture applying even a very small magnetic field rearranges 
the spins so that the spin rotation plane is perpendicu- 
lar to the field direction, in violation of Dzyaloshinskii's 
requirement. Instead of a 2-nd order Dzyaloshinskii tan- 
sition one gets a lst-order spin-flip transition from the 
incommensurate phase directly into the papramagnetic 
state. Although this does not. happen in compounds like 
RbMnBr 3 EH and CsFeCl 3 E], the phase behavior there 
is seriously complicated by the quantum effects and frus- 
tration in the triangular spin latticeEil 

In the following paragraphs we shall demonstrate 
that the experimental data on the CI transition in 
Ba2CuGe207 agrees with our solution to Dzyaloshin- 
skii's model at the quantitative level. Ba2CuGe2C>7 
thus appears to be the first system that exhibits a 
Dzyaloshinskii-type transition in its original form. 

C. Estimates for the critical field 

The appeal of the theory presented above is that it al- 
lows for an exact solution. Its major limitation of course 
is that it is based on a classical, rather than quantum spin 
model. Nevertheless, our theoretical predictions seem 



to be in excellent agreement with what is experimen- 
tally observed in Ba2CuGe207, even at the quantitative 
level. To begin with, the model gives the correct value 
for the critical field H c . The in-plane exchange parame- 
ter J = 0.48 meV was previously determined by measur- 
ing the spin-wave dispersion spectrumjlj. The exchange 
energy per bond is then J = 2JS 2 « 0.24 meV. The g- 
values were measured in ESR experiments: g a = 2.044 
and Qc = 2.474.E3 By substituting C(0) = 0.027 into 
Eq. (E4j) we immediately obtain H c = 3.3 T, that should 
be compared to the experimental value H c ss 2.1 T. Con- 
sidering that for the theoretical estimate we used pre- 
dictions for T — and ignored quantum effects, a 30% 
consistency is quite acceptable. 

D. The soliton lattice 

The most intriguing prediction of our theory for CI 
transition in Ba2CuGe207 is that the spin structure close 
to the phase transition is no longer an ideal spiral, but 
rather should be viewed as a lattice of solitons, i.e., do- 
main walls separating regions of Neel-like spin arrange- 
ment. The soliton lattice is a distinguishing character- 
istic of all CI systems. Several comprehensive reviews 
on the subject (past, among them papers by BakEH and 
Pokrovsky et al.E3 The basic physical mechanism is quite 
simple. When the potential term is sufficiently large, but 
still smaller than the critical value, it is favorable for the 
system to have large "commensurate" regions. The elas- 
tic energy from the "incommensurate" term in the Hamil- 
tonian is partially released by forming domain walls or 
phase slips that separate "commensurate" domains. 

From the experimental point of view the soliton lat- 
tice concept has three important consequences. The 
first is that the incommensurability parameter £ is field- 
dependent, and C — * continuously as H — ► H c (Fig. |^). 
Moreover, the transition itself is very unusual: continu- 
ous, yet logarithmically steep at H c . The experimental 
((H) data collected for the well-aligned domain (Fig. |^, 
open circles) can be nicely fit by the theoretical curve 
shown in Fig. ||, treating ((0) and H c as adjustable 
parameters. A very good agreement is obtained with 
C(0) = 0.027(1) and H c = 2.13(2) T at T = 2.4 K (solid 
line in Fig. |). All the way up to ((H)/((0) ~ 2/3 the 
theoretical curve follows the experimetal points closely. 
In the very proximity of H c however, a deviation is appar- 
ent. This is to be expected: the soliton lattice becomes 
infinitely soft at the transition point and pinning to struc- 
tural defects leads to a saturation of ((H) at H — ► H c 
and effective broadening of magnetic Bragg reflections 
(Fig. ||). The almost-first-order transition is very fragile. 
A small misalignment of the field produces what we see 
in the misaligned domain as a first-order transition with 
a mixed-phase region and a substantially different form 
of C(^0 m the vicinity of H c . 

Inherently related to the unusual logarithmic phase 
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transition is an anomaly in the magnetic susceptibility 
at H c . We have used the theoretical expression ( ^5|) to 
fit the experimental x(-ff ) shown in Fig. ^. Since we can- 
not explicitly take into account the quantum and thermal 
spin fluctuations, we have treated x\\i X-Lj and H c as in- 
dependent fitting parameters. In addition, we allowed 
for a linear term vH in x(H) that is present in both 
the H||a and H||c data. This term empirically accounts 
for intrinsic non-linearities in magnetization curves in 
the quantum (quasi) 2-dimensional AF Heisenber model 
(Rcf. |34|, esp. Fig. 5), that in our case effectively mod- 
ify the local susceptibilities x± an d X\\ & t high fields. 
The fit is shown in a solid line in Fig. g. The values 
of the fitting parameters are: x\\ — 0.89 x 10 -5 emu/g, 
XX = 3.43 x 1(T 5 emu/g, v = 1.72 x 1CT 10 emu/g, and 
H c = 1.88 T. This value for H c obtained from magneti- 
zation at T = 2 K is slightly lower than H c = 2.13(2) T 
obtained from neutron diffraction at T — 2.4K. 

Finally, an important feature of the magnetic structure 
at finite fields is that it no longer is an ideal sinusoidal 
spiral. As discussed above, this distortion is character- 
ized by higher-order Bragg harmonics, that should be 
observable in neutron diffraction experiments. Although 
we have spent some time looking for the 3rd order mag- 
netic satellite during experiment II, at T — 2.4 K and for 
several values of applied field, we were so far unable to 
find it. The reason is probably the lack of intensity. As 
can be seen from Fig. ^ the relative intensity of the 3rd 
harmonic is very small except in the immediate proxim- 
ity of H c , where the 1st satellite itself is broadened and 
weakened. We would like to emphasize here that the form 
of ((H) and the existence of higher-order harmonics are 
in essence one and the same effect. Since theory agrees 
so well with experiment as far as ((H) is concerned, we 
are confident in that satellites are present and will be 
observed in future experimental efforts. 

E. System dimensionality 

We finally comment on our measurements of critical 
exponents. In our previous work we have shown that 
there are several hints, including the temperature depen- 
dence of magnetic susceptibility, to that Ba2CuGe2C>7 
should primarily be considered a 2-dimensional antiferro- 
magnet. This is further confirmed by the measured order- 
parameter critical exponent /3 = 0.184 for the PM^ com- 
mensurate transition, i.e., substantially smaller than in 
standard 3-D models, where f3 > 0.3 for Heisenberg, XY 
and even Ising systems. One could naively expect that 
at high external fields the critical index would be smaller 
than at H = 0, since the symmetry of the corresponding 
Hamiltonian is smaller. Exactly the opposite is observed 
experimentally. The only suggestion that we can make at 
this point is that Ba2CuGe207 falls into a completely dif- 
ferent universality class than conventional magnets. Con- 
sidering the nature of the spiral phase and D-M interac- 



tions, we may— be dealing with one of the chiral univer- 
sality classes£3 To obtain further insight into the crit- 
ical properties of Ba 2 CuGe207 further field-dependent 
measurements of the of (5 and other critical indexes are 
required. 

No matter what the observed critical exponents are, 
the long-range magnetic ordering is still a purely 3- 
dimensional phenomenon. For a CI occurring at fi- 
nite temperature dimensionality is known to play a key 
role. In a purely ID system at non-zero temperature the 
ground state will be destroyed, since the energy required 
to create a soliton is finite while the entropy gain propor- 
tional to InL is infinite in thermodynamic limit. Solitons 
will be spontaneously created at any temperature and no 
sharp transition will occur. In two dimensions, the effect 
of thermal (or quantum) fluctuations is more subtle. For 
a general case Pokrovskii, Talapov and Bak have demon- 
strated that the effective soliton-soliton interaction is al- 
tered by fluctuations.^ The short range exponential re- 
pulsion exp(— Ar) is replaced by a long range term 1/r 2 . 
This modifies the behavior close to the transition point, 
making the phase transition a usual second order-type 
with ((H) ~ (H c — H) 1 / 2 . Only in 3 dimensions the re- 
sults that we have derived for the ground state should 
remain valid at finite temperatures. The transition in 
this case is of almost first order, with logarithmic "cor- 
rections" near the transition point making it continuous. 

F. Ideas for future experiments 

Ba2CuGe207 is a very interesting system, yet it is rela- 
tively easy to investigate experimentally. The Heisenberg 
exchange constants are small, so the magnon dispersion 
relations could be measured in the entire Brillouin zone. 
D-M interactions are relatively strong and the incom- 
mensurability parameter £(0) is sufficiently large to be 
easily measurable. H c 2 T is also readily accessible 
in most types of experiment, even those that require the 
use of diffraction-adapted horizontal-field magnets. The 
work on Ba2CuGe207 is far from being completed. The 
first priority is to find the higher-order harmonics of the 
incommensurate magnetic peaks and study the temper- 
ature dependence of their intensities. One should also 
look for satellites around ferromagnetic zone-centers to 
see if there is a weak-ferromagnet distortion of the spiral 
structure. The effect of magnetic field applied in the (a, b) 
plane is worth investigating. The critical behavior is not 
fully understood and accurate field-dependent measure- 
ments of all critical indexes are highly desirable. Finally, 
it would be interesting to look at dynamical propeties 
of the soliton lattice. Near the critical field the soliton- 
soliton interaction is extremely weak, which gives rise to 
spin waves with very low velocity. 
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VI. CONCLUSION 

We have observed a rare type of CI transition in an 
antifcrromagnetic insulator with Dzyaloshinskii-Moriya 
interactions. The transition is driven exclusively by the 
changing strength of the commensurate potential. The 
latter is directly controlled in an experiment by varying 
the magnetic field. A transition of this kind was envi- 
sioned over three decades ago by Dzyaloshinskii, and we 
have now found that Ba2CuGe207 exhibits it in its origi- 
nal form. In addition, we have extended Dzyaloshinskii's 
theoretical treatment to derive exact parametric equa- 
tions for the field dependence of magnetization and in- 
commensurability vector. 
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FIG. 1. a) Magnetic structure of Ba 2 CuGe 2 7 (Ref. |l5|). 
b) D-M interactions in the Cu-planes of Ba2CuGe207. The 
(l,-l,0)-component of the Dzyaloshinskii vector D (solid ar- 
rows) is the same for all oriented Cu-Cu bond (dashed 
arrow) along the (1,1,0) direction. The ^-component is 
sign-alternating. 



FIG. 2. Insert: Magnetization of single crystal 
Ba2CuGe2 07 measured as a function of magnetic field ap- 
plied along the c and a crystallographic axes. Main panel: 
Field-dependence of magnetic susceptibility x — dM/dH de- 
duced from magnetization curves shown in the insert. The 
solid line represents a fit to the data, as described in the text. 



FIG. 3. a-c) Elastic scans measured along the (1 + C, C,0) 
direction at T = 1.4 K for several values of magnetic field 
applied along the c axis of the crystal, d) Field-dependence 
of the intensity in the (1, 0, 0) Neel peak (solid symbols) and 
the (1 + C,C>0) satellite (open circles). 
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FIG. 4. Measured field dependence of the magnetic propa- 
gation vector £ in the incommensurate phase of Ba2CuGe2 07. 
The data shown by solid circles correspond to a domain mis- 
aligned by fa 1° with respect to the applied field, and were 
collected on a field-cooled sample at T = 1.4 K. Open cir- 
cles are data points measured at T = 2.4 K in an almost 
perfectly aligned domain, obtained by zero-field-cooling the 
sample from T = 5 K. The solid line is a fit with a theoretical 
curve described in the Theory section. 

FIG. 5. Measured magnetic phase diagram of 
Ba2CuGe2 07. The field is applied along the c-axis of the 
crystal. 

FIG. 6. Elastic scans measured along the (1 + f , — £, 0) di- 
rection in Ba2CuGe2 07 at T = 2.4 K for several values of 
magnetic field applied along the c axis of the crystal. For this 
domain the (1,0,0) component does not appear before the 
critical field is reached. 

FIG. 7. Measured Bragg intensity of the (1, 0, 0) (open cir- 
cles, H = 3.66T, T N = 3.5 K) and (1 + CC.0) (solid circles, 
H — 0,Tjv = 3.2 K) magnetic Bragg reflections. The solid 
lines lines are power-law fits to the data. The dashed line is 
an aid to the discussion in the text. 

FIG. 8. Theoretical field-dependence of the magnetic prop- 
agation vector £, obtained from an exact solution of a classical 
1-dimensional antiferromagnet with Dzyaloshinskii-Moriya 
interactions in the continuous limit. Insert: schematic rep- 
resentation of the spiral phase (H = 0), soliton lattice 
(0 < H < H c ) and the spin-flop commensurate state 
(H > H c ). 

FIG. 9. Theoretically predicted Bragg intensity of the 3rd 
Fourier harmonic of the soliton lattice plotted as a function 

ofC(tf)/C(o). 
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